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Abstract — Price of anarchy, the performance ratio, which could 
characterize the loss of efficiency of the distributed supply chain 
management compared with the integrated supply chain man- 
agement is discussed by utilizing newsvendor problem in single 
period which is well-known. In particular, some of remarkable 
distributed policies are handled, the performance ratios in each 
case which have been investigated in the previous works are 
analyzed theoretically and the tighter upper bound of price of 
anarchy and the lower bound are presented. Furthermore our 
approach is developed based on a generalized framework and a 
geometric interpretation of price of anarchy is appeared via the 
literature of convex optimization. 

Index Terms — newsvendor problem, price of anarchy, convex 
optimization, inequalities, geometric interpretation, autoregres- 
sion model with x square noise 

I. Introduction 

MEASURING the efficiency of supply chains plays 
an important role in operations management because 
there exit many complicated factors (various contracts, order 
policies and physical structures) which could influence the 
decision making process. In practice, several formulae have 
been reported in most of the earlier studies [31, 13], E), lfT4ll . 
|fT9l . II20II . In history, bullwhip effect, which could compare the 
variance of orders with that of demands, has been widely used 
to evaluate the efficiency of supply chain management in the 
stationary market |j5]. Although it is straightforward to forecast 
the optimal order in the following term via finite instance 
of demand and order observed according to the manner of 
bullwhip effect, if the market behaves in equilibrium, it is too 
hard to determine the optimal order of goods in the typical 
case with respect to the given supply chain. While, as another 
approach, recently numerical experiments in Cachon (2004) 
indicate that the relative efficiency of a two-stage decentralized 
supply chain could be as low as 70-90% for either push or 
pull configurations with price-only contract policy t4j . Perakis 
and Roels (2007) mathematically extends Cachon's works into 
several different supply chain configurations such as push 
or pull inventory positioning, two or more stages, serial or 
assembly systems, single or multiple competing suppliers, and 
single or multiple retailers lfl4l . By introducing the concept 
of price of anarchy into supply chain, which can measure 
the ratio of the performance of the centralized system to the 
worst performance of the decentralized system (Koutsoupias 
and Papadimitriou 1999 and Papadimitriou 2001), they found 
that even in a two-stage supply chain the loss of efficiency 



might be more than 42% under the same conditions of price- 
only contract, and pointed that a pull inventory configuration 
generally outperforms a push configuration ifTTl . ifTJl . Il20l . 

The contribution of Perakis and Roels (2007) is one of the 
most pioneer investigations in quantifying the performance of 
supply chain. However, since their argument is hard to be guar- 
anteed mathematically, rigorously and sufficiently anywhere, 
we need to improve several points in their logic for more 
practical use. In this paper, tighter upper bound of price of 
anarchy is theoretically proved and a lower bound of price of 
anarchy is presented firstly due to more accurate treatment. 
Furthermore a geometric interpretation of price of anarchy is 
also demonstrated which can give a clear illustration of the 
loss of efficiency of the given distributed supply chain. We 
only consider the both bounds of price of anarchy in two 
echelon case because of the most fundamental case, however 
it turns out that our approach would be simply applied in more 
complicated case. 

This remainder of the present paper is organized as follows; 
in the next section, the newsvendor problem is presented 
for simplicity of our discussion, and the derived results in 
the previous work lfT4l is introduced briefly. Furthermore an 
improved upper bound and a lower bound of the loss of 
efficiency are explained. Section Hill addresses a generalization 
of newsvendor problem and assesses the optimal inventory 
levels and the performance ratios in each mode. In Section 
II VI a geometric interpretation of price of anarchy is intuitively 
provided and we confirm that the performance ratio and the 
analytical procedure handled here are possible to be one of 
the most unbeatable frameworks. The final section is devoted 
to a summary and future work. 

II. Model setting 

Concerned with distribution of goods and information shar- 
ing, supply chain management is one of the most vital interests 
in the cross-disciplinary fields. As one of the most pivotal 
topics, we discuss here the manner of the inventory man- 
agement theoretically, in particular, (1) how to determine the 
optimal inventory level with respect to the given supply chain 
management and (2) how to assess the loss of efficiency of 
suboptimal policy. For simplicity of our argument, we restrict 
the model which is well-defined and is introduced below. 
Therefore applying our demonstration, one could indeed im- 
prove and develop this approach so as to resolve more practical 
case. 
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A. Newsvendor problem 

Newsvendor problem is modeled as follows; if the firm 
prepares the inventory of goods Q and the order in the market 
is ^ , his profit in the single period is expected as follows; 



7r(0 : = -cQ +pmin(Q,^), 



(1) 



where the purchasing cost and the selling price describe c 
and p, respectively, furthermore min((5,^) denotes the lesser 
value of Q and ^. Note that c < p is needed in nature since 
the firm won't prefer to stock and buy the goods in the case 
c > p (briefly the profit 7r(^) < at c > p in other words), 
and notice that Q and ^ are assumed as nonnegative and real 
numbers without the loss of generality. Here the opportunity 
loss is not handled and the inventory space is large enough as a 
matter of convenience, however it turns out that our approach 
would be simply applied in the case with opportunity loss 
and the upper bound of inventory level. Generally speaking, 
it is too hard to estimate that the demand in each single 
period is fixed. Hence, let us propose that the demand ^ is 
stochastically drawn from the given density function /(^) with 
the cumulative probability. Now, the expected aggregate profit 
is represented as follows; 

/.Q _ 

n:= / dU{S.)Ai)^-c^Q+P d^FiO, (2) 



where F(^) dxf{x) describes the cumulative proba- 

bihty whose stochastic variable is greater than or equal to ^ 
and n is a concave function of the inventory quantity strictly 
(show appendix lA-Bb . Note that 11 = at Q = is required 
for any distribution of demand by definition and it implies that 
no stock is no benefit. Moreover as trivial, the supremum of 
the expected entire benefit 11 is indeed greater than or equal 
to zero because of the previous notice. 

The production planner's purpose in general is to maximize 
his expected gross benefit 11 by adjusting the inventory level 
Q in typical situation. However, nowadays the logistics, the 
distribution of goods in supply chain behaves like bloodstream 
in the human society, the necessity not only of the integrated 
inventory management but also of the distributed inventory 
management has been recognized. In the given distributed 
configurations, the problem who makes to store the inventory 
and/or who needs to decide the wholesale price in supply chain 
is one of the most vital issues, furthermore it is also important 
how one examines the loss of efficiency of the distributed 
management in some of remarkable configurations compared 
to the integrated management. 

B. Centralized supply chain 

Let us review here the problem how the inventory level 
Q is derived in order to maximize the expected entire profit 
of the given integrated (or centralized) supply chain. From 
eq. Q, the unique optimal solution which can be desired is 
intuitively derived as follows; Qc := F (r) where r := c/p 
and F (y)(= x) represents the inverse function of F{x){~ 
y). Since the second term in eq. (|2]) is a concave function of 
Q (show appendix lA-BI ). it turns out that the unique optimal 
solution is strictly determined. 



C. The profits in two echelon case 

The counterpart of the centralized case in supply chain 
management, the distributed (or decentralized) supply chain 
management is explained here. In the distributed case the 
expected profit could be divided into two distinguished profit 
functions as follows; 11 = 11^ + 11^ in push serial supply 
chain and 11 = S^' + in pull serial supply chain. Firstly, 
:= {w-c)Q and n^* := -luQ+p d^F {£_) indicate the 
manufacturer's entire profit and the retailer's whole benefit, 
respectively, in the case that the retailer makes to stock the 
inventory in push serial supply chain. While "EJ^ := —cQ + 
w d(F{^) and (p - w) J^^ d(F {^) describe the 

manufacturer's aggregate profit and the retailer's total benefit, 
respectively, in the situation that the manufacturer makes to 
store the inventory in pull serial supply chain. In each case, 
the leader decides the wholesale price w (note that c < w < p 
is required in practice, because these profits are satisfied with 
positivity) in order to maximize the leader's expected whole 
profit, while the follower should choose selfishly the optimal 
inventory level Q by employing the optimization problem 
of the follower's expected benefit with respect to the given 
wholesale price, that is, it is the scenario of Stackelberg 
leadership game ITSl . In the previous investigations, the 
analysis of the class of increasing generalized failure rate 
distribution, where one ensures that the optimization problem 
in the decentralized configurations could possess the well- 
defined optimal solution, has been reported comparatively well 
with regard to several concrete distributed configurations in 
two echelon case as follows; 

(a) The manufacturer is the decision maker in push serial 
supply chain. Qd is determined by the following 
equation; F (Qd) (1 - g{Qd)) = r where g{Q) := 
'^I^Q^ is utilized. In this paper g{Q) is assumed as 
a nondecreasing function of Q and < g{Q) < 1 
because it is guaranteed that the optimal solution of 
the follower's optimization problem is unique. Thus 
g{Q) is termed as increasing generalized failure rate. 
In addition, the desirable wholesale price is yielded 
as u; = pF{Qd). 

(b) The retailer is the decision maker in push serial 
supply chain. Qd is consistent with the inventory 
level in the integrated supply chain management, 
because the inventory is stored at the downstream 
site and it is to be expected that the wholesale price 
is equal to the purchasing cost. 

(c) The manufacturer is the decision maker in pull 
serial supply chain. Qd is also consistent with the 
inventory level in integrated supply chain, because 
the inventory is stored at the leader's site and it is 
to be desired that the wholesale price corresponds to 
the selling price. 

(d) The retailer is the decision maker in pull serial 
supply chain. Qd is decided by the following equa- 
tion; F{Qd){l + l{Qd)y^ = r where 1{Q) := 
—2^^^ So ^^F{£,) is employed and the optimal 

wholesale price is derived as w = c/F{Qd)- If g{Q) 
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is increasing generalized failure rate, then 1{Q) is 
nondecreasing function of Q strictly ID, lfT4l . 
In practice, the production planner should choose the most 
appropriate distributed inventory management in some re- 
markable configurations by some means. For that reason, our 
purpose is here to examine the following measurement so as 
to assess the loss of efficiency of each configuration. Price 
of anarchy, the performance ratio which can characterize the 
loss of efficiency of the expected whole benefit of the given 
decentralized inventory management compared with that of the 
given centralized inventory management and which is well- 
known, is denoted as follows; 



PoA 



-cQ,+p d^FiO 



(3) 



cQd+P d^FiO 



where eq. (|2|i is utilized and it turns out that PoA is greater 
than or equal to unit in general (show appendix I A- Ab . 

D. Rigorous results derived in the previous works 

With respect to the ensemble of increasing generalized 
failure rate, the rigorous results of price of anarchy were 
presented in the previous work as follows llT4l : (a) PoA < 
(1 - k)-^ - (1 - (b) PoA = 1, (c) PoA = 1 and (d) 
PoA < + where fc := g{Qd) and / -.^liQd) 

are utilized (the indices are mentioned in the previous subsec- 
tion). Indeed, although the previous work explained that the 
upper bounds in (a) and (d) are comparatively tight, it is hard 
to assess the loss of efficiency in precision and in good faith, 
therefore we discuss more conscientiously and obtain a tighter 
upper bound compared with the derived one in Perakis and 
Roels (2007), and a tighter lower bound with respect to the 
class of increasing generalized failure rate exactly. 



E. An upper and lower bound of cumulative probability 

From the definition of increasing generalized failure rate, 
logF(^) = — dy^Y^ is analytically yielded where F{0) = 
1. Moreover with regard to the ensemble of increasing general- 
ized failure rate, max (F (QJ , F (Qd) {QdY <F{0< 
max(F{Qc) ,F{Qd)iQd)''C') is obtained where Qd < 
£, < Qc and k ~ g {Qd) and s := g {Qc) are employed. 
According to the ratio of the inventory level in the integrated 
case to the one in the distributed case, a := Qc/Qd, we 
obtain an upper bound and a lower bound of the integration 
as follows; 



< 



C{a,s) < 

i QdF{Qd) 
\ C(a,k) 



d£,F{0 



i-fc 



(1 
(1 



k)' 

k)- 



k < a 

i < a < (1 



k)- 



where C{a, t) QdF {Qd) 

Notice that a < (1 - fc) 
condition, F{Qc) = 



(1 - k)a + 



(4) 



is used. 



is not satisfied with the boundary 
and C{a,t) is a nonincreasing 



and convex function of t for any a. Since the integration 
in eq. (|4|i is evaluated more accurately compared with 
their insufficient discussion in ||T4| . the comparative 
tight upper bound of price of anarchy and the lower 
bound are expected fortunately. Furthermore, since the 
cumulative probability is a nonincreasing function, the 
upper bound in this interval Qd < ^ is not possible to 
exceed the probability at Q = in nature, the inequality 

F{Qd) < F{0 < minfl,F(Qd)(Qd)''r'') is derived, 
and an upper bound and a lower bound of the integration, 

QdF {Qd) < / d^F{0 



< QdF{Qd) 



{Qd, 

1 - k 



(5) 



is obtained. 

F. Both bounds of price of anarchy; the manufacturer is the 
leader in push serial supply chain 

According to the argument in the previous work lfT4l . by 
definition, one can replace price of anarchy as follows; PoA = 



1 



Q di{F(i)~r) 



upper bound 



2 — = — . Thus employing eq. Q and eq. an 




-(l-fc)-i (l-fc)-*< 



" — (1 — fc) ^ otherwise 



and a lower bound 

s{l-k)^-i-s kri-\l-k)^~i ~k 
PoA > - 



(6) 



1 - s 



1 - k 



, , k - kri^'^il - k) 
^+ 1-k 



(7) 



are simply yielded. It turns out that the numerator of the first 
term of the upper bound in (1 — k)~~ < a < (1 — k)^^ 
in eq. (|6ll is a nondecreasing and concave function of a. 
The supremum of the right hand side in eq. ^ was already 
presented in |[l4l, while the comparative tight upper bound is 
to be desired in(l — <Q;<(1 — fc)^^. Fig. [T] shows 

that in the limit of fc for any s > fc, both bounds at 
Q ^ (1 — fc)~~ are close to unit and if s = fc, the lower 



bound 1 



2-k 



1— rfc ' 



is greater than or equal 



to unit at a > (1 — fc)^'^. While Fig. |2] indicates the behavior 
at fc = 0.20 and r — 0.40 and it is found that both bounds are 
greater than unit at any ratio a. 



G. Both bounds of price of anarchy; the retailer is the leader 
in pull serial supply chain 

In this case, compared with the previous subsection, fc and 
s are rewritten as 1 — fc = (1 + l)^^ and 1 — s = (1 + 1)~^, 
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Price of anardiv 




upper bound 



lower bound 



derived bounds are compared with each other in the case that 
the manufacturer is the decision maker in push serial supply 
chain and in the case that the retailer is the decision maker 
in pull serial supply chain, respectively. In conclusion, it turns 
out that our approach is valid in this model. 



3.5 



Fig. 1. The ratio a v.s. botli bounds of price of anarcliy at fc = 0.01, 
^(Qd) = 0.5 and s ~ 1.0 

Price of auarcliv 



1.8 - 




Upper Ijoimd 



lower bound 



2 2.C 
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Q,;, the retailer is the leader hi puh serial supply chain. 
Q,/. the manufacturer is the leader in push serial sup])ly chain 




Fig. 3. Tile ratio r = c/p v.s. Qc and 



Pok 



Fig. 2. Tlie ratio a v.s. both bounds of price of anarchy at k 
FiQd) = 0.5 and s ~ 1.0 



0.20, 



respectively, where / = I (Qd) and t > I, then an upper bound 
and a lower bound; 




PoA < 



PoA > 



(1 + 0'+^ -(1 + {l + l)'+T<a 



Fig. 4. The ratio r v.s. price of anai'chy and the bounds. In push serial supply 
chain, the manufacturer is the decision maker. 



I 



— (l + I) otherwise 



t{l + l)t -t + l-lrT{l + l)T 



{1 + l)-^ + I - IrHl + l)f 
are also yielded. 

H. Example: Nonnegative order drawn from normal Gaussian 
distribution 

Let us confirm the effectiveness of our approach with a 
novel toy model. For simplicity, it is assumed that the demands 
are independently and identically distributed according to most 
of the previous works. Here the density function of demand, 
/(^) is satisfied with — |=e~^ for ^ > and otherwise. 



1.5 



The upper bound derived in the previous work 



The u])per bound 



PoA 



The lower bound 




One can easily validate that g{Q) of this model is increasing Fig. 5. The ratio r v.s. price of anarchy and the bounds. In pull serial supply 

generalized failure rate. Thus and in the two cases of '^^ain, the retailer is the decision maker. 
PoA 7^ 1, are illustrated in Fig. [3] Furthermore, as shown in 
Fig- Hand in Fig. |5] that the numerical results of PoA and the 
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/. Example: Autoregression model with x square noise 

As the real-world data analysis, it is natural that the de- 
mands are correlated with one another, for instance, a market 
trend, rather than that they are independently and identically 
distributed ifTOll . ifTTl . ifTSl . However, in most of the earlier 
studies, for simplicity, the demands were independently and 
identically distributed and this assumption is obviously not 
practical since real-world date are usually somewhat biased 
and correlated across the instance. In general, the demand 
might consist of the trend effect which implies correlation and 
an uncertain factor which connotes noise (c.f. bullwhip effect). 
Therefore as the first step, let us handle a novel toy model 
in order to confirm whether or not these bounds of price of 
anarchy which have been presented in this paper are influenced 
by the correlation in the market. Well, the demand at discrete 
time T (denoted by ^t) is generated from the following 
autoregression model; ^t+i = /^Ct + c^Xt' where (3 indicates 
a dumping coefficient in < /3 < 1 and x^, the random 
variable at time T, is drawn from x square distribution with 
one degree of freedom (in addition, for simplicity each noise 
is assumed to be independently and identically distributed 
and implies a noise intensity). Although there is not an 
aggressive premise, because the demand is guaranteed to be 
always nonnegative, this toy model is accepted here. Now, 



0.01 




(3 = 0.9 and cr^ = 100.0 in the case that the manufacturer is 
the leader in push serial supply chain and in the case that the 
retailer is the leader in pull serial supply chain, respectively. 
Likewise, it turns out that our procedure is validly supported 
in this correlated model. 



2500 



Qfj, the retailer is the leader in pull configuratiou. 



manufacturer is the leader in push configiuation. 




(1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
r 

Fig. 7. The ratio v.s. the inventory levels at /3 = 0.9 and = 100.0. 
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:).0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Fig. 8. The ratio r v.s. price of anarchy and the derived bounds. In push 
serial supply chain, the manufacturer is the decision maker. 



Fig. 6. It is well-known that the density function /(§) is analytically yielded 
via the cumulant generating function, the logarithm function of characteristic 
function, log rf5/(?)e'«« = -i J2T=i '°S " 2461^2/3'=) 0. How- 
ever we estimate here the density function utilizing the histogram numerically 
evaluated by 0.5 X 10^ demands which are consisted of the sequences which 
are randomly chosen from the times series, which is stable, and we fit the 
logarithm of the density function derived numerically into the fifth-degree 
polynomial function of log 5, which is supported by minimization of the the 
leave-one-out cross validation error |18| . /(^) at /3 = 0.9 and cr^ = 100.0 
is illustrated in this figure. 



given (3 and cr^, and the density function of demand /(^) 
is stable, one can evaluate the loss of efficiency, PoA with 
respect to the ratio r = c/p, furthermore, we can confirm 
whether or not both bounds are valid. The density function of 
demand, /(^) at /3 = 0.9 and = 100.0 is illustrated in Fig. 
|6l Qc and Qd in the two cases are depicted in Fig. [7] As shown 
in Fig. [S] and in Fig. |9] that the numerical result of PoA and 
the upper and lower bounds are compared with one another at 



PoA 




Fig. 9. The ratio r v.s. price of anarchy and the derived bounds. In pull 
seiial supply chain, the retailer is the decision maker 
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III. Generalized Framework 

We have discussed price of anarchy, the performance ratio 
which can characterize the loss of efficiency of the given 
distributed supply chain management compared with the in- 
tegrated management policy, based on newsvendor problem 
which is well-known. In particular, some of remarkable dis- 
tributed policies are handled, the performance ratios in each 
case, for instance, who makes to store the inventory and/or 
who needs to decide the wholesale price, are analyzed the- 
oretically, and the tighter upper bound of price of anarchy 
and the lower bound are presented. However, two points 
should be noted here that these results are not restricted 
to newsvendor problem (hereafter it is termed as original 
newsvendor problem) and price of anarchy is possible to be 
regarded as one of the most unbeatable feature quantity with 
respect to a broad class of optimization problems. 

In this section, let us propose a generalization of the 
approach treated in the previous section (hereafter generalized 
newsvendor problem is handled intuitively, however it is also 
demonstrated in the literature of supply chain management 
as a matter of convenience). According to our argument, 
one could replace the profit function in eq. ([T]i as 7r(^) = 
—cQ + pm{^,Q), i.e. m{^,Q) = min(^,(3); it turns out 
that m(^, Q) can describe a more complicated policy of price 
contracts, which could sufficiently characterize a market trend, 
for instance, the order of hard-to-find items in practice. Now, 
for simplicity of the discussion, we prepare the expected 
aggregate profit. 



n 



= -cQ+pM{Q), 



(8) 



where M{Q) := dUiOH^.Q) and X{Q) 
are defined. In order to extend the method, we should explain 
J^{Q) and X{Q) simply. Firstly m{£,,Q) is here 
an arbitrary function which is satisfied with the following 



conditions. 



d^MjQ) _ dX{Q) 



dQ 



< and M{0) = 0. Because 



it is required for the concavity of the expected total profit 
function (one would not presume that Q) is a concave 
function of Q in general, i.e. m(^, Q) = min(^, Q)) and 
M{0) = is necessary so as to connote that no stock is 
no benefit. Furthermore from X{Q) and M{0) = 0, M{Q) 
is possible to be replaced as A4{Q) ~ d^'^{^) (it is not 
required to disclose the density function of demand in gen- 
eralized newsvendor problem), where X{£^) can be regarded 
as F{£^) in the previous section. Further it implies that X{Q) 
is assumed as positive because the optimal inventory level is 
nonnegative. Lastly, c < pX{0) is needed in nature (briefly 
n < for any inventory level is satisfied at c > pX{Q)). 



(iv) := {p — w)A4{Q), where the wholesale price w is 
satisfied with c < w < pX{0) in push serial supply chain 
and c/X{0) < w < p in pull serial supply chain, respectively. 
According to the previous argument, price of anarchy is also 
denoted as follows; 



PoA 



pMiQc) 



-cQd+pMiQd)' 



(9) 



In the follows, we explain how the optimal inventory level 
and the desirable wholesale price of each decentralized case 
are determined so as to measure the loss of efficiency. 



B. Manufacturer is the leader in push serial supply chain 

It is well-known that PoA in this situation of original 
newsvendor problem is not always equal to unit. In practice 
one can derive the optimal solution iteratively via the saddle 
point equation as follows; 



Qd = X-\e) 



r + giQd)X{Qd) 



(10) 



where e :— w/p and a novel function, g{Q) := 
~Q-^logX{Q) are already used. By definition, g{Q) is a 
nonnegative function of Q in general. Furthermore in order 
to determine the unique solution of the decentralized manage- 
ment, we assume that g{Q) is increasing generalized failure 
rate (note that if g{Q) is termed as increasing generalized 

failure rate, then g{Q) is satisfied with < g{Q) < 1 and 
dg{Q) > n\ 
dQ - 

In addition, although this iteration connotes the recursive 
procedure in order to resolve the optimal inventory level 
systematically in the distributed system, using eq. ( fTOb . one can 
represent also the following relation as X{Qd) (1 — g{Qd)) = 
r, where the desirable wholesale price is ui = pX{Qd) by the 
definition of e. 



C. Retailer is the leader in push serial supply chain 

In this case, PoA is equal to unit because the derivative of 
the leader's expected profit 11^ with respect to the wholesale 
price w is nonpositive, the optimal wholesale price is desirable 
to be consistent with the purchasing cost in c < w < pX{Q), 
therefore it is possible to be regarded as the integrated system 
as the follower's benefit is zero, that is, Qd is equivalent to 

Qc 



A. Centralized case and profit fimctions 

The unique optimal solution in the integrated case of 
generalized newsvendor problem is denoted as follows; Qc ~ 
X~^(7'), where r = c/p (notice that < px(o) — 1) and 
X^^{y){~ x) represents the inverse function of X{x){= y). 
Next, the expected profit function 11 in eq. ([8]l is possi- 
ble to be divided into two distinguished parts in push and 
pull configurations as follows; (i) II'^ := {w — c)Q, (ii) 
n« -wQ+pM{Q), (iii) := -cQ + wM{Q) and 



D. Manufacturer is the leader in pull serial supply chain 

Fortunately, PoA coincides with the value derived in origi- 
nal newsvendor problem because the derivative of the leader's 
expected profit S'^ with respect to the wholesale price w is 
nonnegative, the optimal wholesale price is desirable to be 
equal to the selling price in c/X{0) < w < p, therefore 
it is possible to be regarded as the integrated system as the 
follower's benefit is zero, that is, Qd is consistent with Qc- 



7 



E. Retailer is the leader in pull serial supply chain 

It is sure that PoA in the last case of original newsvendor 
problem is not less than unit. With respect to generalized 
newsvendor problem, let us evaluate the optimal solution 
sequentially utilizing the steepest descent method as follows; 



evaluated as follows; 



1 



KQd) 

XiQdV 



(11) 



where 5 := clw and a novel function liQ) := —-^-t^^-ttM: 
are employed. Under the definition, 1{Q) is a nonnegative 
function of the inventory level in nature. Moreover, so as 
to derive the unique solution of the distributed management, 
we require that 1{Q) is a nondecreasing function of Q as 
the sufficient condition (if g{Q) is assumed as increasing 
generalized failure rate, then 1{Q) is a nondecreasing function 
of Q, show appendix lA-Db . 

Additionally, although this iteration derived here implies the 
algorithmic procedure so as to assess the optimal inventory 
level in the decentralized system, applying eq. ( fTTT i. the 
following relation, X{Q(i) (1 + l{Qd))^^ = '% is obtained 
exactly where the desirable wholesale price is w — c/ X{Qd)- 



F. Both bounds of price of anarchy 

From the above discussion, it turns out that PoA in the two 
distributed cases that the retailer is the leader in push serial 
supply chain and the manufacturer is the leader in pull serial 
supply chain, respectively, namely the inventory is stocked at 
the leader's site, is similar to the derived results of original 
newsvendor problem. Herein the others are discussed. 

Firstly we explain the case that the manufacturer is 
the decision maker in push serial supply chain. Utiliz- 
ing the definition of increasing generalized failure rate, 
logA'(^) = logvf(O) — dy^^ is assessed. Thus with 
respect to the class of increasing generalized failure rate 
distribution, max(A'(Q,),-Y(Qd)(Qd)"e"') < ^(0 < 
max{X{Q^),X{Qd){Qd)'T'') is derived where Qd < ^ < 
Qc and k = g{Qd) and s = g{Qc) are represented, respec- 
tively. According to the ratio a = Qc/Qd, one can analyze an 
upper bound and a lower bound of the integration as follows; 



JQd 

- I Cia,k) 



(1 - k)-i < a 
(1 - k)-i <a<il~k)- 



where C{a, t) 



QdX{Qa 



(1 - k)a + 



is employed. Moreover QdX{Qd) 



< 



1-t 



< 



QdXiQd) 



1 



k 

1-k 



1 



is calculated since 



X{Qd) < < mm{X{0),X{Qd){&r'') is obtained 

in ^ < Qd- Therefore an upper bound and a lower bound 
of price of anarchy in generalized newsvendor problem are 



PoA < 



PoA > 



(l-fc)--^ -(l-fc)-i (l-fc)--^<a 



I, fc(i-fc) 
s(l-fc)i- 



- ( 1 - fc) - 1 otherwise 
(l-k) 



kri-^l-k^-i -k 



1 - s 



1 



, , k-kr-!^-\i-ky-- 

^ + l-k — 

where we replace the ratio r as the rescaled ratio of the 
purchasing cost to the selling price, r := ^;y{o) ^ Pjl]- It 
turns out that the derived bounds in generalized newsvendor 
problem are as well as the ones in original newsvendor 
problem. Lastly, both bounds are also simply derived in the 
case that the retailer is the leader in pull serial supply chain. 

G. Example: A toy model 

The previous argument has indicated only if one vahdates 
both bounds of the performance ratio in resolving generaUzed 
newsvendor problem, we need not to restrict to the literature of 
supply chain management. As a matter of course, with respect 
to the ensemble of increasing generalized failure rate without 
the context of operations management, one also needs to vin- 
dicate the improved and developed bounds. Hence we would 
apply M.{Q) := tanh((5) for simplicity of the discussion 
because X{Q) = 1 - tanh^(Q), g{Q) = 2gtanh(Q) and 
1{Q) = 2sinh^((5) are derived briefly and analytically. 

Firstly, in the case that the manufacturer is the leader in push 
serial supply chain (as a matter of convenience we address so), 
it is comparatively easy to execute the algorithm based on eq. 
([Tol l, then Qc and Qd are illustrated in Fig. [TO] and PoA and 
the derived bounds are indicated in Fig. [TT] Conclusionally, it 
turns out that our approach is valid in this case. While, in the 
case that the retailer is the leader in pull serial supply chain, it 
is also to perform the iteration based on eq. ( fTTl i. then Qc and 
Qd are shown in Fig. [10] and PoA and the derived bounds are 
presented in Fig. [T2] Similarly it turns out that our procedure 
is supported validly in this case. 

IV. Geometric interpretation 

In the previous section, we have examined price of anarchy 
with respect to generalized newsvendor problem. Furthermore 
here, in order to comprehend the measure in depth |6|. let us 
provide a novel geometric interpretation with respect to price 
of anarchy via the property of the convexity of —M.{Q) ||2l, 

ESI. 

A. Geometric interpretation; the integrated supply chain 

First, one can divide two distinguished functions with 
respect to the expected aggregate profit function, 11 = — cQ + 
pM (Q) as follows; j/i (Q) = cQ + H and 2/2 (Q) = pM {Q) 
where from yi(Q) — y-ziQ), that is, if there exists intersection 
point, then 11 = —cQ + pAi (Q) is derived. As shown in 
Fig. [T3]that 11 implies the intercept of yi{Q). Or as another 
representation, the linear function of Q which has slope c 
and passes through an intersection point (Q* ,pM (Q*)) is 
represented as follows; yl (Q) = c{Q — Q*) + pM{Q*). 
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dyi (Q) dv2 (Q) 




0.1 0.2 0.3 0.4 0..5 0.6 0.7 0.8 0.9 



Fig. 10. The ratio r = c/p v.s. Qc and Qd- 




Fig. 11. The ratio r v.s. price of anarchy and the bounds. In push serial 
supply chain, the manufacturer is the decision maker. 



PoA 




The upper bound 



The lower bound 




PoA 



01 02 03 04 05 06 07 08 09 



Fig. 12. The ratio r v.s. price of anarchy and the bounds. In pull serial 
supply chain, the retailer is the decision maker. 



Indeed the intercept of yl{Q) describes also H. Show Fig. [T3] 
in order to maximize the intercept of yi{Q), both functions 
should intersect at one point (denoted by {Qc,pA4 (Qc))) 
at least. Since the derivatives of both functions of Q, i.e. 



dyiiQ) 



c and 



dy2{Q) _ 



dQ 



dQ 



at Q 



Qc, namely 

pM (Qc) , 



(12) 



is possible to be sufficiently satisfied. Thus we can resolve the 
unique optimal solution of the integrated inventory manage- 
ment. 




{Q'.pM{Q-]] 

\i{Q)=pM(Q) 



Fig. 13. The maximum of intercept of linear function implies the extremum 
of the expected total profit function. 



One point should be worthy to be noticed here. X (Q), the 
derivative of Ai (Q) of Q, is not always necessary for the 
continuous function of Q (however by definition, A4{Q) is 
satisfied with the continuous function because of the concav- 
ity). For example, the given function. 



pM (0) 



log (1 + 0) 

v{Q - Qv) 



< < 

l0g(l + g.„) Qy<Q 



dQ 



= pX (Q), are yielded easily. 



is defined with constant Qy and v where v < jqj^ is required 
because X{Q) is a nonincreasing function of Q. Then with 
respect to the slope of yi{Q) in ?; < c < , it turns out 
that the intersection point is {Qv,pM {Qv)), however, this 
optimal solution is not satisfied with eq. ( fT2] i indeed. That is, 
we need to comprehend that eq. (fT2l) describes the sufficient 
condition but not the necessary. Without the loss of generality, 
so as to prevent also us from misleading in practice, we should 
confirm the behaviors of both functions yi{Q) = cQ + 11 and 
2/2 (Q) = pM{Q) being supported by the picture such as Fig. 

m 

B. Geometric interpretation; the manufacturer is the leader in 
push serial supply chain 

According to the explanation in the previous subsection, 
with respect to 11^ = —wQ + pA4{Q), in order to solve 
the optimization problem of the follower, the two functions, 
yi{Q) ■= wQ + and y2{Q) ■= pM{Q) are denoted. 
As shown in Fig. [14] that the supremum of the intercept of 
yi{Q) (the optimal solution is {Qd,pM {Qd})) describes the 
desirable whole profit and the expected total profit in the 
centralized case, 11 (Qc) is greater than maxll'^ because of 
c < w. Furthermore since the leader's aggregate profit is 
represented as II'^'' = (w — c)Qd, Fig. [14] illustrates that 
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- n{Q,) = -cQ, + pM (Q.) 




(Qd.pM(Qi)) 


"A 


<l^/y2{Q)=pM{Q) 




- max = ~wQ,, + (Q,,) 




Q 



Fig. 14. A geometric interpretation of price of ana rchy; the manufacturer 
is the leader in push serial supply chain. Fig. 1151 illustrates how Qd is 
determined. 




Q 



Fig. 15. Qd is satisfied with pXiQd) (1 - g(Qd)) = P^(Qc). 



n iQc) > {w-c)Qd+i-wQd + pM (Qd)) =: n (Qa), where 
Il{Qd) connotes the expected profit of the decentralized case. 
Conclusionally, it turns out that PoA is greater than unit in 
nature. 

C. Geometric interpretation; the retailer is the leader in pull 
serial supply chain 

In the other case, so as to analyze the optimization problem 
of the follower, with respect to S'^ — —cQ + Wj\4{Q), the 
two functions, yi{Q) := cQ + E^ and y2{Q) ■= wM{Q), are 
defined. As shown in Fig.[T6]that the maximum of the intercept 
of yi{Q) (the optimal solution is {Qd,pM (Qd))) represents 
the desirable entire profit and the expected aggregate profit in 
the integrated case 11 {Qc) is greater than maxS^ because of 
w < p. Moreover if the leader's whole profit is = (p — 
w)MiQd), Fig. [I6] depicts, n{Q^) > -cQd + wM (Qd) + 
{p — w)A4 (Qd) = Il{Qd)- In conclusion, it turns out that 
PoA is greater than unit in general. 

V. Conclusions 

We discussed price of anarchy, the performance ratio, 
which could characterize the loss of efficiency of the dis- 



Il{Q, ) = ~''Q, +pM (Q,,) 




Q 



Fig. 16. A geometric interpretation of p rice of anarchy; the retailer is the 
leader in pull serial supply chain. Fig. [IT] illustrates how Qd is determined. 




Q 



Fig. 17. Qd is satisfied with pX(Qd) (1 + KQd))"^ = P^iQc). 



tributed supply chain management compared with the inte- 
grated supply chain management via newsvendor problem and 
the generalization instead of bullwhip effect. In particular, 
the performance ratios in some of remarkable decentralized 
supply chain managements are analyzed theoretically and 
numerically. Furthermore, with respect to the ensemble of 
increasing generalized failure rate, which one can ensure that 
the optimization problem of the follower could possess the 
well-defined solution; (a) the upper bound which has been 
investigated in the previous work ifTTI is improved in this 
paper utilizing the more accurate evaluation of the integration 
and (b) the lower bound is derived in the same manner, in 
the case that the manufacturer could control the wholesale 
price in push serial supply chain and in the case that the 
retailer could adjust the wholesale price in pull serial supply 
chain, namely the two cases that the follower makes to stock 
the inventory. Moreover the framework handled in section 
has been developed and deepened for generalized newsvendor 
problem, and we indicate that the loss of efficiency is measured 
as well as original newsvendor problem. Hence price of 
anarchy is possible to be one of the most unbeatable feature 
quantity with respect to the convex optimization involving 
Stackelberg leadership game lfT6l . While our approach is 
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supported validly in some examples which are satisfied with 
increasing generalized failure rate. Lastly without the loss of 
generality, a geometric interpretation of price of anarchy has 
been provided concretely. 

The investigations of geometric interpretation of price of 
anarchy in multiechelon case, and of the other ensembles 
which are guaranteed that the optimization problem of the 
follower can possess the well-defined solution are promising 
topics for future works. 
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Appendix A 
Preliminaries 

A. Global optimal and local optimal 

In this appendix, let us introduce the relationship between 
global optimal and local optimal. Well, we assume that 
X and y are convex sets and employ as x ^ X and 
y E y, respectively. Furthermore two real-valued functions 
f{x,y) and g{x,y) are bounded from above in the region 
(a;, y) G X ^y. Then a novel function is denoted as follows; 
F{x,y) := .f{x,y) + g{x,y) where this function is also 
satisfied with bounded above and let {x*,y*) be an extremal 
solution of F{x, y) in the given finite region. While [x** ^y**) 
indicates an extremal solution of f{x,y), that is, one part of 
F{x, y), then, F{x*,y*) > F{x**,y**) is obtained in general. 
Thus price of anarchy is greater than or equal to unit by 
definition. 



B. The concavity of d^F{£,) 

Because F(^) is a nonincreasing function of ^ firstly, 

Jo d^FiO < Jo" d(FiO+FiQo){Q-Qo) is held for any Q 
and Qo in general, therefore, for any Q, Q' and I G [0, 1], the 

concavity, I d^O + (1 - /o'^' d^O < Io^T{0^^ 
satisfied when Qo = + (1 — l)Q' is rewritten. Additionally, 
if Q — Qd and Qo = Qc, it is also proved that PoA is 
greater than or equal to unit from d^F{^) < jQ"d(F{^) + 
F{Qo){Q~Qo). 



C. Young's inequality 

Because F(^) is a nonincreasing function of ^, Qip + 
J^dyF~^{y) > Jq dCF{^) > QF{Q) is held without the 
loss of generality for any f > Q. The more left inequality 
is termed as Young's inequality, iff F{Q) — if gives the 
equality. Therefore the more right inequality is obtained from 
Jo d(F{i) - QF{Q) = dyp- {y) > at F(Q) = ip 

without the property of increasing generalized failure rate Q, 
0. 



D. A proof of nondecreasing function 1{Q) 

If g{Q) is increasing generalized failure rate, then 1{Q) is 



strictly satisfied with 



> 0. Because one can prepare 



a novel function firstly, j{Q) := a'(q) /o^ d^X{^), where 
= 1 + 1{Q). Here we allow to rewrite 1{Q) as 

m 

to Q 



j(Q)g{Q) ^ pj-om the derivative of log/(Q) with respect 



1 dl{Q) _ 1 dg{Q) 



9{Q) - 1 - 



jiQ) 



1{Q) dQ g{Q) dQ Q 

is calculated where Q, j{Q), g{Q) and 1{Q) are nonnegative 
by definition. Hence as proof by contradiction, g{Q) < 1 — 

-^H^ is assumed. Then the derivative of 1 — rfL^ is derived 

j(Q) ](Q) 

to be negative exactly. Thus g{Q) > > 1 — is yielded 



1 



where limg^o g{Q) = limg^o , x 
result is inconsistent with the assumption g{Q) < 1 
namely g{Q) > I - 

proved that 1{Q) is a nondecreasing function of Q 



1(0) 

= 0, however this 



Q 



j[Q) ' 

is held in nature. Therefore it is 



E. Magnitude relation of price of anarchy 

From Young's inequality and the discussion of the previous 
appendix, 1{Q) > g{Q) and l{Q)g{Q) - 1{Q) + g{Q) > are 
derived, respectively. Therefore the relationship between the 
derivative of the inventory level Qd^puii in pull serial supply 
chain with respect to the rate r and the derivative of the 
inventory level Qd.push in push serial supply chain with respect 

dOaj^ < Strictly. Thus 



to the rate is obtained as 



, dr 



Qd.puu = -Irdr^%^ > -J'dr^^f^ = Qd,push is 
held. Moreover because 11 is a nondecreasing function of Q 
in Q < Qc, PoA in push configuration is greater than or equal 
to PoA in pull configuration in general. 

F. Price of anarchy in the fixed order case 

For any Q G (0, Qc], the sufficient and necessary condition 
of the equality 1 - g{Q) = (1 + l{Q)y^ is M{Q) oc Q. 
In original newsvendor problem, it implies that F{£^) ~ 1 for 
^ < Qo and 0, otherwise, namely the order is fixed at ^ = Qq, 
then Qc = Qd = Qo is held in each decentralized supply 
chain. In conclusion, if the order is constant, the performances 
of the four cases discussed in this paper are consistent with 
one another. 

Appendix B 

n serial supply chain management 

Our approach based on generalized newsvendor problem is 
simply to be extended N serial supply chain management (the 
case of iV = 2 is already mentioned). The optimal inventory 
level of each distributed management is devoted as follows: 

(a) The manufacturer is the decision maker in push sup- 
ply chain. Qd is derived from the following equation; 

(l + Q-^) "^(Q) = ^ where r = c/p and 
roughly speaking, Qd is probably to be satisfied with 
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the condition X{Q) (1 — g{Q))^~^ > r. Here one 
point should be noteworthy. The optimal inventory 
level Qd is not always satisfied with the equality, 
X{Q) (1 - g{Q)f'^ = r, because = is 

not always held for each integer in 1 < ?i < 
Moreover it is hardly desirable that the comparative 
tight both bounds of price of anarchy in this con- 
figuration are derived by the optimal inventory level 
which is satisfied with the above inequality. 

(b) The retailer is the decision maker in push supply 
chain. Since the inventory is stored at the leader's 
site, Qd ~ Qc is desirable. 

(c) The manufacturer is the decision maker in pull 
supply chain. Nevertheless to say, as the goods is 
stocked at the upstream site, Qd = Qc is expected. 

(d) The retailer is the decision maker in pull supply 
chain. Qd is derived as the solution of the following 

relation + -J- 
relation, | i H ;f(Q) dQ ) 



where 



roughly speaking, Qd is possible to be satisfied with 
the condition X{Q) (1 + 1{Q)) 



-{N-l) 



> r. 



Appendix C 
Multiple materials and multiple items 

We could develop our approach in the case of the in- 
ventory management of S multiple materials and / multiple 
items briefly. Let c := {ci,C2, - ■ ■ ,cs}'^ S and Q := 
{Qi, Q2, ■ ■ • 1 Qs}^ G R-'^ t>6 the purchasing costs and the in- 
ventory levels of the materials, respectively. Furthermore p := 
{Pi,P2,---,Pif G and M := {Mi,M2,---,Mif £ 
represent the selling prices and the order levels of the 
items, respectively. The entry of the order levels is assumed 
to be strictly a convex function of Q. The expected whole 
profit is defined as follows; 11 = ~c^Q + M., where the 
notation T denotes the matrix transpose. First, the optimal 
inventory levels of the integrated supply chain Qc is held with 
the following equation. 



E 



dQs 



Next, in the case that the manufacturer is the leader in push 
serial supply chain, the inventory levels of the decentralized 
configuration Qd is satisfied as follows; 



I 

Cs = 



dQs 



dQsdQt 



. Lastly, 



Therefore the wholesale prices w := {tui, W2, • 
is denoted as follows; Ws = ^i=iPfj. (^^) 
in the case that the retailer is the decision maker in pull 
serial supply chain, the inventory levels of the distributed 
configurations Qd is satisfied as follows; 



s 

= E' 

s=l 



dQs 



d^Qs 



i/=i 



dMndM^ 



where the desirable wholesale prices w is held as follows; 

Qd 



w, 



Three points should be noted here. Firstly, several specific 
problems which are possible to be handled via our approach 
have been already investigated |[T], fS), lfT2l . Especially, only if 
the solution with respect to the given optimization problem is 
definite in several types, utilizing such as Lagrange's multiplier 
method, we could resolve in the same way. Next, in the 
case of the supply chain management in multiperiod (i.e. 
S = I), let Cs, Ps, Qs and be the purchasing cost, the 
selling price, the inventory level and the demand at term s, 
respectively. Further the demands ^ := {^1, ^2, ■ ■ ■ , Cs}^ are 
distributed with the given multivariate density function /(^) 
and Ms /g°° d£,f{^)ms{^, Q) indicates the order at term s 
where ms(^, Q) describes that the order at term s is depended 
on the market trend strongly. Likewise, it turns out that one 
can deal with multiperiod case based on our approach (c.f. 
bullwhip effect). Lastly, in the case of N serial supply chain, 
the optimal inventory levels Qd in the two cases are satisfied 
with the following equations; 



Cs = Y^Pt. 

M=l 

S 

= E 



7V-1 



1=1 



Tr n [5u^^.+M,^- 



dQs 



respectively, where t :— {ti,t2, ■ ■ ■ ,tN-^i} and j7 := 
{1^1, 1^2, • • • , '^N-i} are used, 5ab indicates Kronecker's delta 
which is the entry of the unit matrix and the notation Tr^- and 
Tri7 denote the summation over all possible states of t and of 
ly, respectively. 
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